Coarsening of Precipitation Patterns in a Moving Reaction-Diffusion Front 
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Precipitation patterns emerging in a 2D moving front are investigated on the example of NaOH 
diffusing into a gel containing AICI3. The time evolution of the precipitate Al(OH)3 can be observed 
since the precipitate redissolves in the excess outer electrolyte NaOH and thus it exists only in a 
narrow, optically accessible region of the reaction front. The patterns display selfsimilar coarsening 
with a characteristic length £ increasing with time as £(£) ~ \ft. A theory based on Cahn-Hilliard 
phase- separation dynamics including dissolution is shown to yield agreement with the experiments. 
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Understanding precipitation patterns formed in a mov- 
ing reaction front is important both from the basic aspect 
of extending our knowledge of phase separation dynamics 
and from the point of view of technological applications. 
Indeed, the formation of precipitation structures, well lo- 
calized in space and time, underlies the notion of the so 
called " bottom up" designs p], [H, 0, 0] where one creates 
structures directly in the bulk. A natural way of real- 
izing such design is to have a reaction-diffusion process 
and control the dynamics of the reaction zone i.e. control 
both the position of the front and the rate of the creation 
of the reaction product which yields the precipitate pro- 
vided its concentration exceeds a threshold value. 

Recently, there have been several attempts at control- 
ling reaction zones [3|, [4|, |5|, |6[ with the most promising 
results emerging from the use of the ionic nature of the 
reagents and realizing control through time-dependent 
electric currents [7j . Although these experiments demon- 
strate that one can create one dimensional structures re- 
producibly at the scale of ~ 500/i, the downsizing of the 
patterns to submicron scales raises several problems re- 
lated to the front. One of them is the width of the front 
which obviously restricts the downscaling in the direc- 
tion of the motion of the front. The second one is related 
to the inhomogeneities within the front which may lead 
to unwanted precipitation structures in the plane per- 
pendicular to the motion of the front. The scale of these 
structures clearly limits downscaling in the transverse di- 
rection. 

The width of reaction fronts has been extensively 
studied theoretically @, 0] as well as experimentally 
OH EH 03 and one has some ideas about the control 
parameters in this case. Very little is known, however, 
about the transverse patterns in a moving reaction zone, 
though bulk coarsening has been studied in connection 
with the so called gradient-free precipitation experiments 

iQ. 

Our aim here is to initiate experimental and theoreti- 
cal studies of the coarsening dynamics of the transverse 
patterns in reaction zones. This task is made feasible ex- 



perimentally by overcoming the transparency problem in 
a way suggested by earlier studies of Liesegang type phe- 
nomena [iBl, EB, E3, El, [10, [2l|| . Namely, the reaction- 
diffusion process is set in a nearly transparent gel and, 
furthermore, appropriately chosen electrolytes are used 
so that the reaction product undergoes redissolution in 
the excess of the outer electrolyte. As a result, precipitate 
exists only in a narrow region restricted to the reaction 
zone and its time evolution can be followed in detail. 

We observed the patterns in the moving reaction zone 
in an experimental setup detailed below. The visual ob- 
servations suggested that the system displayed a selfsim- 
ilar coarsening and this was quantified through the time 
dependence of the characteristic lengthscale of the 
pattern, with the main experimental result being that 
~ y/t. Theoretically, the effective dimensionality 
of the coarsening system is not entirely obvious, and 
there are several candidates for driving the coarsening 
process. We studied this problem by generalizing the 
Cahn-Hilliard theory of precipitation to include sources 
and sinks coming from the emergence of the reaction 
product in the reaction zone, and from the redissolution 
of the precipitate in the excess outer electrolyte, respec- 
tively. The numerical solutions of the equations in three 
dimension are in agreement with the experimentally ob- 
served £ (t) ~ \ft. This suggests the natural picture that 
the sources and sinks are relevant perturbations on the 
particle conservation, and we observe a curvature driven 
late-stage coarsening in a model with nonconserved order 
parameter [22|, [23[ . 

In the experiments, a 1 w/w% agarose (Reanal) gel 
was prepared with height of 0.6 — 0.7 cm in a Petri dish. 
After the gelation process took place (« 2— 3 h), the inner 
electrolyte (AICI3) was poured on top of the agarose gel 
to obtain a given concentration in the gel (0.48 mol/L - 
0.56 mol/L). After 64 hours the inner electrolyte solution 
was removed from top of the gel and replaced by the outer 
electrolyte (NaOH of fixed concentration 2.5 mol/L). 

A white precipitation layer formed immediately at the 
gel interface, and this layer started to move into the gel. 
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The evolution of precipitation pattern in the moving layer 
was recorded in reflected light using a EOS-20D camera 
connected to a computer. Typical recording period was 
£ — 25 minutes and we analyzed the processes in this time 
window supposing that the light intensity is proportional 
to the precipitation concentration [241 ]. 
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FIG. 1: Time evolution of the precipitation pattern in the 
reaction zone for the samples with [NaOH] = 2.5M (outer 
electrolyte) and [A1C1 3 ] = 0.52M (inner electrolyte). The 
front moves perpendicularly to the plane of the picture and 
the pictures were taken at £1 — 180s (a), £2 = 480s (6), and 
£ 3 = 960s (c) after the initiation of the reaction. The length 
of the scalebars is 1cm. 

As shown in FigJIJ the patterns have a random appear- 
ance and they coarsen with time. The coarsening displays 
selfsimilarity as indicated by Fig|2] where the pictures 
have been magnified by l/y/i with £ being the elapsed 
time from the initiation of the process. The selfsimilar- 
ity can be quantified by measuring the time evolution of 
the characteristic length £(£) of the precipitation struc- 
tures. We found £ 2 through calculating the structure 
factor of the grayscale values of the pattern and evaluat- 
ing the average of the wavenumber squared £ 2 ~ l/(k 2 ). 
The results for three different inner electrolyte concentra- 
tions are displayed in Fig J3] and we can see a well defined 
diffusive growth £ 2 ~ 2D±t regime. 




FIG. 2: Visual demonstration of the scaling in the evolution 
seen in Fig[T] Parts of the panels (a) and (b) are magnified by 
factors (£s/£i) 1//2 and (£3^2) respectively. The scalebars 
are the same length of 1cm. 

The transverse diffusion coefficient D± can be esti- 
mated from FigO and we find D± = £ 2 /2£ « 5 • 10~ 10 
m 2 s _1 . It is remarkable that D± is an order of magni- 
tude smaller than the diffusion coefficients of the small 
hydrated ions. 

Our understanding of the observed phenomena is as 
follows. The hydroxide ion (outer electrolyte) diffuses 
into the gel, and the white precipitation layer at the gel 



FIG. 3: Square of the characteristic length £ 2 of the pre- 
cipitation patterns plotted as a function of time £ for three 
inner electrolyte concentrations and at fixed outer electrolyte 
concentration (00 = 2.5M). We calculated £ 2 from the 2nd 
moment of the structure factor of the patterns. 

interface is produced by the reaction with the inner elec- 
trolyte (aluminum ions) 

Al 3+ (aq) + 30H-(aq) Al(OH) 3 (s) . 

Next, this layer redissolves due to the complex formation 
of aluminum hydroxide in the excess of hydroxide ions 
producing a soluble aluminum complex 

Al(OH) 3 (s) + OH-(aq) -> [Al(OH) 4 ]"(aq) . 

The combination of the precipitation and of the complex- 
formation processes results in a thin precipitation layer 
moving through the gel. It should be emphasized, how- 
ever, that the motion is solely the transport of chemical 
species and not of the precipitate. 

Within the above picture, the coarsening in the reac- 
tion zone emerges from a complicated interplay of the 
dynamics of the reaction zone, of the aggregation onto 
the already present precipitate, and of the redissolution 
process. As a first step in understanding the interplay 
of these processes, we approached the problem through 
coarse grained (mean field) equations by using a model 
[25( | that has been developed in a series of papers during 
the last decade 0,13. 

The three processes we include are as follows. First, 
the reaction of the electrolytes A + B — > C yields the re- 
action product C. Since the process takes place in a gel, 
no convection is present and it can be modeled as a sim- 
ple reaction-diffusion process. This reaction provides the 
source for the precipitation which is modeled as a phase 
separation of Cs described by the Cahn-Hilliard equation 
with a source term. Finally, the Cs redissolution in the 
excess As (complex formation; A + C — > Complex) ap- 
pears as a sink term both in the Cahn-Hilliard equation 
and in the reaction-diffusion equation for the As. 

The picture can be further simplified if we assume that 
A + B — > C is an irreversible process for totally dissoci- 
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FIG. 4: Simulation results for the spatial distribution of C integrated along the z direction. The parameters used were: 
bo/ao = 0.18, c h /a = 0.2, c z /a = 0, D = 10" 9 mV\ ki = 8.6 M - ^" 1 , k 2 = 0.13 M'V 1 , A = 0.83 x 10" 9 mV\ a = 10" 6 
m 2 , the grid spacing and the time step were 0.01 mm and 0.003 s, respectively. The snapshots are taken at 75s, 600s and 
1200s. 



ated electrolytes A and and that the case of monova- 
lent ions with equal diffusion coefficients are considered. 
Then the equations take the following form: 



dta — DAa — k\ab - 
d t b = DAb-hab 
d t c = -XA(5f/5c) 



koac 



k\ab — k^ac + \fc 



(1) 
(2) 

crj. (3) 



Here D is the diffusion coefficient of the ions, while k\ 
and &2 are the rates of reaction and of complex formation, 
respectively. We take k\ to be large resulting in a reac- 
tion zone of negligible width (note that this assumption 
is compatible with the typical reactions used in experi- 
mental setups that result in Liesegang structures). The 
thermal fluctuations in c are described by a noise term 
^fcr\ which conserves the total number of C particles [26[ . 
The free energy (/) underlying the thermodynamics of 
the phase separation ([3]) is assumed to have minima at 
some low (q) and high (ch) concentrations, and it is as- 
sumed to be of a Ginzburg-Landau form. Its functional 
derivative most compactly given in terms of a shifted and 
rescaled concentration m = (2c — Ch — q) / (ch — q) 



Ch ~ ci Sf _ 6f_ 
2 5c 5m 



m — m + a Am . 



(4) 



Finally, a and A in eq.Q are setting the spatial- and the 
time scales. They can be chosen to reproduce the correct 
time and lengthscales in experiments [27| . 

The initial conditions to the above equations are set 
according to the experiments. The outer and inner 
electrolytes are homogeneously distributed in the lower 
(z < 0) and upper (z > 0) halfspaces with concentra- 
tions ao and bo <C ao, while the initial concentration of 
C is Co = everywhere. For solutions in finite rectangu- 
lar boxes, we applied no-flux boundary conditions at the 
borders along the z axis while periodic boundary condi- 
tions were used at the rest of the border planes. 

Equations (Q][3|) together with (jH were discretized on 
a uniform 3D grid of size 144 x 144 x 704, and finite- 
size effects were checked on grids of size 160x160x568 



(no effects were observed) and 96x96x1584 (effects were 
seen: the final characteristic length was increased roughly 
by factor of two). 

The resulting ordinary differential equations were in- 
tegrated in time by Euler method (fast simulations were 
made feasible by the parallel programming possibilities 
of video cards [28|] ) . The results from the simulations are 
summarized in Figs. [4] and [51 

Fig. [4] suggests the existence of two stages in the evo- 
lution of the system. There is an initial period when the 
unstable, homogeneous density of Cs produced by the 
front begins to evolve towards the equilibrium densities 
(Panel a in Fig. [4j). This stage is driven by the initial 
perturbations (noise) which determines whether the con- 
centration in a given neighborhood grows or diminishes. 
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FIG. 5: Simulation results demonstrating the diffusive growth 
£ 2 rsj t of the characteristic lengthscale £. 

The first stage is finished when the high- or low- 
concentration states are reached in a significant frac- 
tion of the system and an initial phase of fast coarsen- 
ing has also taken already place. The emerging high- 
concentration regions (Panel b in Fig. 0]) form a struc- 
ture which is the base for the next stage, the self-similar 
coarsening. It should be noted that the memory of this 
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structure can be recognized throughout the later stages 
of the coarsening (Panel c in Fig. 

In order to compare with experiments, we calculated 
the structure factor and the characteristic length £(t) by 
processing pictures as in case of the experiments. The 
results for various inner electrolyte concentrations are 
shown in Figj5j As we can see, there is a fast initial stage 
governed by noise and short-scale relaxation to the equi- 
librium concentrations, followed by the late-stage coars- 
ening. The initial large fluctuations in £ are due to the 
fact that, at the initial stage, there is a significant prob- 
ability for two precipitation bands to coexist. The aver- 
aging of concentration in the z-direction yields then an 
apparent pattern with random correlation length. Since 
the Cahn-Hilliard equation describes the long-time, long- 
wavelength properties of the coarsening, the short-time 
fluctuations of £ may be an artificial feature of the re- 
sults. Similarly, the explanation of the presence of an in- 
duction time in the bo = 0.48 M experiments may also be 
outside the scope of the CH description. The late stage 
coarsening, however, should be correctly given. is 
indeed smooth in this regime and, in agreement with the 
experiments, it shows diffusive behavior (£ 2 ~ t). 

We should emphasize that the emergence of a coarsen- 
ing state is by no means obvious. The rate of the arrival 
of A and B particles to the front decreases with time 
(roughly as a§j \]t and bo / y/i) and the front can advance 
only in the presence of surplus As. Thus it is also a con- 
ceivable scenario that, due to the redissolution process 



(A + C — > Complex), there is only small concentration of 
Cs in the front with finite (or perhaps even decreasing) 
correlation length. The selection of a selfsimilar coarsen- 
ing state is the result of a delicate interplay between the 
diffusive advance of the front and the reaction-diffusion- 
aggregation processes within the front. 

In conclusion, we studied the coarsening of precipita- 
tion patterns in a thin moving reaction front, and our 
experiments suggest that the asymptotic dynamics of 
the system may be interpreted as curvature-driven, late- 
stage coarsening in systems with nonconserved order pa- 
rameter. The theoretical approach, based on the Cahn- 
Hilliard equation coupled to reaction-diffusion processes, 
reproduces all relevant findings observed in experiments, 
but the question of why selfsimilar coarsening is selected 
by the dynamics is still to be answered. We believe that, 
in general, the coupling of reaction fronts with phase sep- 
aration processes opens a wide range of possibilities for 
studying new aspects of coarsening and pattern forma- 
tion, and developing this field is of importance for sub- 
microscopic technological design. 
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